Trirefringence in nonlinear metamaterials 
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We study the propagation of electromagnetic waves in the hmit of geometrical optics for a class 
of nearly transparent nonlinear uniaxial metamaterials for which their permittivity tensors present 
a negative principal component. Their permeability are assumed positive and dependent on the 
electric field. We show that light waves experience triple refraction - trirefringence. Additionally to 
the ordinary wave, two extraordinary waves propagate in such media. 
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Introduction. — Electromagnetic waves in a mate- 
rial medium propagate according to Maxwell's equations 
complemented by certain relations linking strengths and 
induced fields - the constitutive relations. Depending on 
the dielectric properties of the medium and also on the 
presence of applied external fields, a variety of optical 
effects can be found. One of wide interest is the birefrin- 
gence [H-Q , occurring when electromagnetic waves prop- 
agate in media exhibiting two distinct refractive indexes 
[l| in a same wave vector direction. The phenomenon 
of triple refraction has been much less investigated. By 
trirefringence in a given wave vector direction, we mean 
the existence of three distinct refractive indexes in that 
direction. In the realm of linear electrodynamics, trire- 
fringence does occur Nevertheless, nothing was stud- 
ied when the dielectric tensors are field dependent (non- 
linear electrodynamics). 

The production of novel optical materials with multire- 
fringent properties has been measured in tailored pho- 
tonic crystals fsj. However, in these media, structural 
details (lattice constants, defects, etc) are imperative and 
hence they can not be described in terms of effective di- 
electric tensors 0]. Developments in photonic band gap 
materials and the so-called metamaterials have enabled 
the discovery of several new phenomena [gl . For instance, 
it was experimentally shown Q that nearly transpar- 
ent isotropic metamaterials allow light propagation when 
both effective dielectric coefficients (permittivity and per- 
meability) are negative. In fact, this phenomenon and 
other unusual properties displayed by isotropic media 
with negative dielectric coefficients were long ago pro- 
posed theoretically Q . It is worth emphasizing these me- 
dia must be dispersive and the negative coefficients are 
obtained for convenient frequency ranges. Some other 
unusual properties exhibited by specific metamaterials 
are negative refractive index [9i] , trapping of light [13, U | , 
perfect lens devices [12| , the electromagnetic cloaking ef- 



fect |13l . [14| and the occurrence of asymmetry for pro 



agation of light in opposite wave vector directions 
Wave propagation in indefinite metamaterials (where not 
all the principal components of the dielectric tensors have 



the same sign) has been also considered [i^ , showing that 
effects already proposed in the context of isotropic meta- 
materials can be obtained and possibly improved. In- 
definite metamaterials can also be used for investigating 
certain aspects of General Relativity [16]. 

In this Letter, we show that nonlinear metamateri- 
als described in terms of effective dielectric tensors, may 
display trirefringence. Analytical expressions describing 
this effect are formally obtained from Maxwell's electro- 
magnetism and a simple theoretical model is numerically 
examined. A possible experimental realization of the me- 
dia expected to display this effect is also addressed. The 
vectorial three dimensional formalism [17| is used. The 
units are set such that c = 1 . 

Wave propagation. — The electrodynamics of a con- 
tinuum medium at rest in the absence of sources is gov- 
erned by the Maxwell field equations 
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taken together with the constitutive relations between 
the fundamental fields E and B, and the induced ones D 
and H, written here as 
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The dielectric coefficients and /i^ are the components 
of the permittivity and the inverse permeability tensors, 
respectively, and they encompass all information about 
the electromagnetic properties of the medium. Further, 
for any vector a we denote its ith component by {i = 
1,2,3). 

The propagation of monochromatic electromagnetic 
waves is here examined within the limit of geometrical op- 
tics [l[ using the method of field disturbances [ll] . This 
method can be summarized as follows [l3|- Let E, de- 
fined by (j){t, x) = 0, be a smooth (differentiable of class 
€"■,71 > 2) hypersurface. The function is understood to 
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be a real- valued smooth function of the coordinates {t, x) 
and regular in a neighborhood U oi T,. The spacetime 
is divided by E into two disjoint regions , for which 
(p{t, x) < 0, and corresponding to 4>{t, x) > 0. The 
discontinuity of an arbitrary function f(t, x) (supposed 
to be a smooth function in the interior of ) on S is a 
smooth function in U, and is given by [isj 



where we defined the traces 



[f{t,x)]^= Jrn [/(P+)-/(p-)], 



{P±} 



(4) 



with P+, P~ and P belonging to U^, U~ and E, re- 
spectively. The electromagnetic fields are supposed to be 
smooth functions in the interior of and U~ and con- 
tinuous across E (</> is now taken as the eikonal jlj] of the 
wave). However they have a discontinuity in their first 
derivatives, such that [isj 



ujbi 
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where and bi are related to the derivatives of the elec- 
tric and magnetic fields on E, and correspond to the com- 
ponents of the polarization of the propagating waves [l9| . 
The quantities uj and qi are the angular frequency and 
the ith component of the wave vector. (Incidentally, we 
note that the negative signs appearing in Eq. ([6]) are 
missing in the corresponding equations in Ref. ^17].) 

For the cases of interest in this work, the permittivity 
of the media under study will be described by real di- 
agonal tensors (losses have been neglected) whose com- 
ponents Eij are dependent only upon the constant fre- 
quency of the wave. We set the magnetic permeability 
of these media to be real functions of the modulus of 
the electric field, such that iJ,ij{\E\) = 6ij/^{\E\), where 
dij = diag (1,1,1). Thus, applying the boundary condi- 



tions stated by Eqs. I^H^ to the field equations ((T][2]) we 
obtain the eigenvalue equation j20|, [2l| 



0, 



(7) 



where the Fresnel matrix Zij is given by 
fi' ^ 1 

Zij = Eij 2 ^ 2 (I'^^ij ^ lilj) (^) 
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in Z2 = / Zii Z 



Z3 = / Zii Zji Zh 



(10) 

As a requirement of the geometrical optics limit, the wave 
fields are considered to be negligible when compared with 
the external fields. Thus, we assume from now on that 
the fields are approximated by their external counter- 
parts Eext and Bext- We set E « E^xt ^ Ex and B « 
Bext = B y, which could be arbitrary functions of space 
and time coordinates. Let us examine the particular case 
of uniaxial media where Sij = diag (ey, ej_, ej_). Us- 
ing Eqs. ([5]) and (ITUl) . straightforward calculations show 
that Eq. ([9]) results in the following algebraic fourth de- 
gree equation for the phase velocity v = u / q oi the prop- 
agating waves. 



dv 



0, 



with 



a — 6e± Ey, 
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c= e^[2e^fqx^ + {s^+e^^)iqy^ + q,^)], (14) 



d= :r£±EBqz, 

6 

e= — [e^qx^ +e^{qy^ +q,^)] 



(15) 
(16) 



We defined qa = {q ■ a), where q = q/q, for any unit 
vector a. Then, q — qxX + qyi) + q^z and q^ = q ■ q — I. 
Notice that when the propagation occurs in the xy-plane 
(spanned by the external electric and magnetic fields), 
the coefficients b and d are null and the generalized Fres- 
nel equation reduces to a quadratic equation in v'^, there- 
fore allowing only birefringence. The same behavior oc- 
curs if /i' = 0. In fact in this situation we recover linear 
electrodynamics. 

Solving Eq. (fTTj) we obtain 



1 



(17) 



with Ai' = {l/\E\){dfi/d\E\) and q^ = q- q. 

Non-trivial solutions for the eigenvalue problem in Eq. 
([7]) can be found if, and only if, det | Zij \= 0, which 
is known as the generalized Fresnel equation. This equa- 
tion gives also the dispersion relations of the media under 
study. Using the covariant formulae for the traces of lin- 
ear operators 0] the Fresnel equation can be cast as 



det I Z,, | = -i(Zi) 



— Zo — — Z'^ — , 
2 3 



= -<yq. ± J {<^qzy + ^ (^^qi + ql + ql , , 



where 



^i' EB 
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The solution Vo does not depend on direction of the 

(9) 

wave propagation and will be called the ordinary wave. 
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whereas v"^ depend on direction of wave propagation and 
will be called extraordinary waves Q . By definition, the 
velocities of the waves are given hy v — vq, where v is 
given by Eqs. ((T7)) and ^TE\\ . In order to achieve more 
simplicity in the following analysis we assume the exter- 
nal fields to be constant and set the wave vector in a 
given direction of the xz-plane, i.e., (jy = 0, Qx = sin9 
and Qz = cos 9. In this notation 9 indicates the angle 
between q and the z direction. 

Trirefringence. — Two distinct solutions for in a 
same given direction q can be obtained from Eq. (|18l) if 

- 1 < — ^ ( ^ tan^ 9 + 1] < 0. (20) 

Thus, taking into account the ordinary wave, Eq. (j20p 
defines a region inside which trirefringence occurs in any 
chosen direction q. Let us examine this effect closer. In 
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FIG. 1: (color online). Normal surfaces 0, 0] of a non- 
linear medium with dielectric coefficients given by Sij — 
diag(— ey, e±, e±) and i-i{E). We set a model where the nu- 
merical values were taken such that ey = 6.69, e± = f.88, 
fi = 0.8 and fi'E B = 4.21. The ordinary wave is represented by 
the circular solid line and the extraordinary waves are represented 
by the dashed and dot-dashed curves. The dashed straight lines 
encompass an angular region in which trirefringence occurs. In 
this case there exist two extraordinary waves denoted by + and 
— . The symbol o refers to the ordinary wave. 

order to guarantee the existence of an ordinary wave we 
must set > 0, otherwise Vo is not real. This is true 
when both coefficients and e± present the same sign, 
which can be positive for usual media or negative for 
left handed materials (in this case a negative refractive 
index occurs Billl). Let us set these coefhcients to 
be positive. Now, in order to satisfy Eq. ([20)1 we set 
— — < 0. Hence, trirefringence occurs in directions 



determined by 

— > tan^ 9>—{l-e, ^a^) . (21) 

fill ^11 

The above discussed phenomenon is displayed in Fig. ([1]) 
where the normal surfaces [H [1] associated with the or- 
dinary and extraordinary waves are depicted for some 
specific values of the quantities appearing in Eq. (j2ip . 
For any given direction encompassed by the angles be- 
tween the two dashed straight lines, defined by Eq. (|2ip . 
there are three distinct solutions: the dashed and dot- 
dashed curves representing the extraordinary waves, and 
the circular solid curve representing the ordinary wave. 
For angles between the dashed and the solid straight lines 
only birefringence occurs, and finally only one refraction 
occurs for directions encompassed by the angles between 
the two solid straight lines. It is also worth noticing 
from Fig. ([T]) that in the sectors where more then one 
refractive index occur, the medium under consideration 
behaves as a positive or negative medium [J, Q , depend- 
ing on sub-sectors and extraordinary waves. 

In geometrical optics the directions of light rays are 
given by the directions of the group velocities u — 
duj/dq = vq + q{dv/dq), which are considered as the 
physical velocities of propagation of the rays [ij. As we 
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FIG. 2: (color online). Ray surfaces [T, '5| of the nonlinear 
medium considered in Fig. [T] The ordinary group velocity is 
represented by the circular solid line and the extraordinary group 
velocities are represented by the dashed and dot-dashed curves. 

see from Eq. ([T5|. the extraordinary waves depend on 
the wave vector. Thus, the directions of the extraordi- 
nary light rays do not in general coincide with the di- 
rections of the extraordinary phase velocities. Taking 
V = {vo,v^ ,v^} in the definition of the group velocity 
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we obtain that the ordinary group velocity is identified 
with the ordinary phase velocity Uo = Vocj- Nevertheless, 
the extraordinary group velocities are given by 



having trirefringence derived in this letter will still hold 



UxX + UzZ, 



with 



■ sin 29 



ME II \e± J 



sni f cos 



+ (1 — r])ave cos 6 + cos^ 9 
cos 9 + (Tvl cos 261 - ^( ^ - i] 

_f M£|| \e±- ) 

+ (1 — T])av^ cos 9 + 2^ cos^ I 



cos V sm 



(22) 



■ (23) 



,(24) 



and where we define r\ = {uj/e\\){de\\/duj). In these equa- 
represents the extraordinary wave veloci- 



tions Vf 



ties in the direction which can be obtained from Fig. 
[T] by taking the distances from the associated points on 
its wave vector surfaces to the origin of the coordinate 
system. For the particular model set in Fig. [T]the cor- 
responding ray surfaces [Ull] are depicted in Fig. [2] as a 
function of the angle Lp between u and the z axis, where 
it was assumed for simplicity that \ri\ « 1. Notice that 
this assumption resides in the realm of nonlinear media. 
In the context of linear electrodynamics, for which — 0, 
Eqs. (I22II24I1 show that I77I << 1 leads to superluminal 
group velocity solutions. Therefore, this assumption is 
not in disagreement with the so-called causality require- 
ment 0. For the ordinary waves we have = 9, as 
anticipated. For the extraordinary waves, the relation- 
ship between these angles can be obtained directly from 
Eqs. ((23ll24| a,s tp = arctan[M^(6l)/M;,(6')]. From Fig. [21 
we see that there is a sector inside which there exist two 
extraordinary rays in any chosen direction of observa- 
tion. For the complementary sector there are none ex- 
traordinary rays. We point out that any chosen direction 
for the rays inside the trirefringent sector is related to 
two different extraordinary wave vectors in Fig [TJ Sim- 
ilarly, for any given wave vector inside the two dashed 
straight lines in Fig. [U there are two different extraor- 
dinary rays associated with it in Fig [2] For instance, 
let us take the wave vector in the direction 9 = 7r/16, 
as indicated by the arrow drawn in Fig. [TJ As is clear, 
in this direction there are three solutions for the phase 
velocity, indicating that trirefringence occurs. These so- 
lutions are represented by the points on the +, — and o 
curves, found in this particular direction, and their mod- 
uli are v'^ = 0.76, v~ = 0.20 and Vo = 0.82, respectively. 
From Eqs. (|22II24|) . the corresponding group velocities 
present magnitudes = 0.83, u~ = 0.55 and Uo = 0.82. 
Their associated directions are indicated by the arrows 
appearing in Fig. [2l where the dashed and dot-dashed 
arrows correspond to and u~ , respectively. 

Final remarks. — If the principal permittivity com- 
ponents are also dependent upon the modulus of the elec- 
tric field, then analogous calculations leading to the equa- 
tions for the phase and group velocities and conditions for 



by identifying e,, ^ e^^{E) + E'^[e^^{E)]' , 



XE). 



Once we work in the regime of high frequencies (re- 
quired by geometrical optics) losses must be assumed to 
be low m . Metamaterials with low- losses are still a mat- 
ter of current investigations [i^] , with some achievements 
already obtained in such regime [25| . 

In what concerns the role played by the nonlinearities 
Eqs. p8)) - p9l) show us that if they exist (i.e. ^' ^ 
0), irrespective of their strengths, then they will lead to 
the presence of two extraordinary waves (hence allowing 
trirefringence). Given the physical quantities appearing 
in Eqs. (111])"(1M1), the thresholds for the strengths of 
the nonlinearities are such that causality is not violated. 
The important quantity related to the nonlinearities is cr, 
given by Eq. ()19|) . which is experimentally controllable, 
thereby making the predictions here derived feasible. 

Layered media 0, HI] seem good candidates for exper- 
imental realizations of our assumptions. Consider a sys- 
tem constituted by a repetition of two thin low-loss lay- 
ers. One of them is composed by a nonmagnetic (/ii = 1) 
and dispersive (ei = ei(w)) medium. The other one is 
composed of a liquid medium whose permeability and 
permittivity are dependent upon the modulus of the re- 
sultant electric field as ^2 = l^/d-S'l) and £2 = l—g{\E\), 
where f(|-E|) and (?(|i?|) are usually much smaller than 



unity 27|. Additionally, if convenient external fields are 



present, then it is always possible to guarantee that our 
constraints are fulfilled. 

Refraction analyzes were disregarded in our work once 
they are just a straightforward extension of the analyzes 
valid for birefringent crystals • The only difference now 
is the existence of two extraordinary waves. 

Closing, trirefringence is a phenomenon allowed to oc- 
cur in some specific nonlinear metamaterials in terms 
of effective dielectric tensors and in principle could be 
formulated both in terms of wave and ray propagation. 
With the present technology in manipulating the dielec- 
tric coefficients in metamaterials ^6|, , we hope that 
the effect here derived can be experimentally tested and 
if verified could lead to applications. 
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